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1. Introduction 

The eponymous Kadomtsev-Petviashvili (KP) equation 

(-4ut + u xxx + 6uu x ) x + 3(3uyy = , (1) 

discovered in 1970, describes the dynamics of small-amplitude, long wavelength, solitary 
waves in two dimensions (xy-plane) [Q]], and arises in the study of water waves, plasma 
and various other areas of physical significance (see e.g., J2l for a review). In equation 
(OQ), u = u(x, y, t) is the wave amplitude, the subscripts denote partial derivatives with 
respect to x,y,t, and = ±1. Throughout this article, equation ([TJ) with = 1, which 
corresponds to the negative-dispersion KP equation (KP II) will be considered, and will 
be referred to as the KP equation. 

The KP equation is a completely integrable system whose underlying mathematical 
properties have been extensively studied during the last four decades. They are well docu- 
mented in several monographs including (but not limited to) SIHEH^. These properties 
include the existence of multi-soliton and periodic solutions and the Lax representation of 
the inverse scattering transform. A major breakthrough in the KP theory occurred in 1981 
when Sato 10] formulated the KP equation in terms of an infinite dimensional Grassmann 
manifold known as the Sato universal Grassmannian. A finite dimensional version of the 
Sato theory corresponding to the real Grassmannian Gr(iV, M) (the set of iV-dimensional 
subspaces of IR A/ ) leads to a simple algebraic construction of a special class of solitary 
wave solutions of the KP equation called the line-soliton solutions. These are real, non- 
singular solutions which decay exponentially in the xy-plane except along certain direc- 
tions. Specifically, such a solution is localized along two distinct sets of rays (referred to 
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as line solitons) as y — > ±00, and form spatial interaction patterns in the finite region of 
the xy-plane. 

The simplest example of a line-soliton is the one-soliton solution of KP given by 

u(x, y,t) = a sech 2 (k • r + ujt + 6q) , (2) 

which is a traveling wave with r := (x,y), amplitude a = \{k 2 — ki) 2 , wave vector 
k := (k x , k y ) = \(k 2 — k\, k 2 — k\) and frequency uj = — kf), where kx, k 2 are 
distinct real parameters with k\ < k 2 . Clearly, the soliton amplitude depends on the wave 
vector k, which together with the frequency u satisfy the soliton dispersion relation 

4uk x = 4k* + 2,k 2 y . 

The solitary wave-form given by (f2]) is localized in the xy-plane along a line which makes 
an angle ifi measured counterclockwise from the y-axis where 



tan^ = k y /k x = k\+k 2 , —^<tp< 



A one-soliton solution is shown in Figure 1(a). Since the one-soliton solution is character- 
ized by two real parameters kx, k 2 , it is convenient to denote this solution simply as the 
[1, 2] -soliton. Note that when kx + k 2 = 0, the solution in Q becomes y-independent and 
reduces to the one-soliton solution of the Korteweg-de Vries (KdV) equation. 

The soliton interactions of the KP equation were originally described via a 2-soliton so- 
lution with a "X"-shape pattern in the plane formed by the intersection of two line solitons 
but with a parallel shift of the two lines at the intersection (the phase shift). This 2-soliton 
solution (see Figure 1(c)) referred to as the O-type soliton ("O" stands for original) and 
its iV-soliton generalization were obtained independently by several authors using integral 
equations [8] and other algebraic methods |@,|25l]. In 1977, Miles fllll pointed out that the 
O-type 2-soliton solution becomes singular if the angle of the intersection is smaller than a 
certain critical value. As the angle approaches the critical value from above, the 2-soliton 
phase shift tends to infinity, and at the critical angle the O-type solution degenerates to 



a "Y-shape" with only three line solitons interacting resonantly (see also II12I1 ). It turns 
out that such Y-shape interacting wave-forms are also exact solutions of the KP equa- 
tion lfl3l. [l4ll . Apart from the ones mentioned thus far, no other soliton solutions of the 
KP equation were known for quite some time until recently when more general types of 
resonant and non-resonant line-soliton solutions were reported in several works including 
1 15l 16l 17l Till . During the last 5 years, considerable progress has been made toward the 



problem of classifying all exact line-soliton solutions of the KP equation 



These studies have revealed a large variety of soliton solutions which were totally over- 
looked in the past. Generically, these solutions of KP consist of two distinct sets of line 
solitons of different amplitudes and along different directions in the xy-plane as y — > ±00. 
Locally, each line soliton denoted as the [i, j] -soliton has the form of a one-soliton solu- 
tion as in (0, parametrized by two distinct real parameters ki < kj, i < j. The soliton 
amplitude is given by c&ij = \{kj — ki) 2 and the soliton angle ipij (counterclockwise from 
the y-axis) satisfies tan^- = ki + kj. 

There are several direct and algebraic methods to construct the KP line-soliton solutions 
which are usually derived from either a Wronskian or from a Gram determinant (Gram- 
mian). In principle any given soliton solution can be represented in either form, although 
this fact has been explicitly shown for only the one-soliton and the O-type soliton solu- 
tions. For more general line-soliton solutions, it is often more difficult to impose regularity 
conditions on the solution when it is represented in one particular form than the other, thus 
rendering one solution generating method less efficient than another. In this note, we dis- 



November 12, 2009 3:17 Applicable Analysis aa 



3 




-100 - 50 50 100 -100 -50 50 100 -100 -50 50 100 



(a) (b) (c) 

Figure 1. Line-soliton solutions of KP: (a) a one-soliton solution with (^1,^2) = (0.5,1) at t = 0, tj> is the soliton 
angle; (b) a Miles's Y-shape soliton solution with (ki,k2,ks) = (—1,0, J) at t = 0; (c) an O-type 2-soliton with 
(fci,...,A;4) = (— 1, — i, 2) at t = 0. In each figure the dominant exponentials in the r-function (see section 2) in 
different regions of the xj/-plane are labeled by the indices in parentheses, and the asymptotic line solitons are labeled by 
the index pair in square brackets. 



cuss both the Wronskian and the Grammian forms of the general line-soliton solutions, 
and establish the equivalence between the two representations in an explicit fashion. That 
is, we show how to derive one form of the solution from the other, and vice-versa. An- 
other purpose of this paper is to describe the interaction properties of the various types of 
line-soliton solutions. In particular, we present a detailed discussion of two distinct types 
of two-soliton solutions both of which form an "X"-shape pattern on the xy-plane but 
interact in a significantly different way. 

The paper is planned as follows. In Section 2, we introduce the r-function of the KP 
equation and present its Wronskian and Grammian forms generating the line-soliton solu- 
tions. Then the equivalence of these two representations of the r-function is established. 
Section 3 is devoted to a brief description of the distinct types of the 2-soliton solutions, 
followed by a detailed discussion of the nonlinear interaction properties of the O-type 
and P-type 2-soliton solutions. We conclude the paper with a brief summary and possible 
significance of the results. 



2. The KP r-function 



The most convenient representation for the line-soliton solutions is via the r-function 
which plays an important role in the mathematical theory of the KP equation 0,0]. The 
solution u(x,y,t) of (Q]) can be expressed as 

u(x,y,t) = 2(lnr) xx , (3) 

where r(x, y, t) is the KP r-function which is defined up to an exponential factor that is 
linear in x, y and t. 



2.1. The Wronskian form of the T-function 

The r-function can be expressed as the Wronskian determinant 0, H3, |25|] 



r{x,y,t) = Wr(/i,...,/jv) 



fi 1 ' f: 



\Jn 



(i) 



e (N-l) 
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where /W denotes the i th partial derivative with respect to x, and where the functions 
{fn}n=i form a set of linearly independent solutions of the linear system 

fy = fxx i ft = fxxx ■ (5) 

A remarkable fact is that the KP equation simply turns into a determinant identity if one 
substitutes (H) into dD and uses the Wronskian form of r(x, y, t) given by (0]) (see e.g., 



la, |23l |25D for details). This implies that any linearly independent set of solutions of the 
linear system (f5]) will give rise to a solution of the KP equation. Hence it is possible to 
generate a large class of solutions in this way. In particular, the line-soliton solutions are 
obtained from the choice 



M 



f n (x,y,t) = 2^ a nm e m , n = 1, 2, . . . , N , (6) 

m=l 

where 9 m (x, y, t) = k m x + k^y + k^t + 0Q m with distinct real parameters: k\ < &2 < 
. . . < ku and real constants {0om}m=l* The coefficients (a nm ) =: A define an N x 
M constant matrix which is of rank N due to the linear independence of the functions 
{fn}n=i- Then © can be expressed as 

r(x,y,t) = ^2 A(mi,...,m N ) exp[ 0(mi, . . . , m N ) ] (k mr - k m J , (7) 

l<mi<— <m N <M l<s<r<N 

by expanding the Wronskian determinant using Binet-Cauchy formula. In above, 
6 (mi, . . . ,m N ) := 9 mi + 6 m2 + ... + 9 mN , and A(mi, . . . ,m N ) is the maximal 
minor, i.e., the determinant of the N x N sub-matrix of A obtained from columns 
1 < mi < • • • < mj^ < M. Note that the N linearly independent rows of the coef- 
ficient matrix A span an iV-dimensional subspace W of IR A/ so that A can be regarded as 
a point of the real Grassmannian Gr(N, M). A different choice of basis for W amounts 
to performing row operations: A — ► GA, G £ GL(N, IR), which changes the r-function 
in {/J simply by a scale factor: r — ► det(G)r but leaves the KP solution u = 2{hxr) xx 
invariant. Consequently the coefficient matrix A can be canonically chosen in the re- 
duced row echelon form (RREF) which has a distinguished set of pivot columns such 
that the restriction of A to this set corresponds to the N x N identity matrix. Given any 
N x M matrix A of rank N, its unique RREF gives a standard coordinate for the cor- 
responding point in Gv(N, M). It is in this sense that the space of all solutions of KP 
generated by the r-function given by ((TJ) can be identified with the real Grassmann mani- 
fold Gr(A r , M). In general, such solutions are not all line-solitons as they can be singular 
where the r-function becomes zero. The regularity of the solutions in the entire xy-plane 
and for all values of t is achieved by imposing a further restriction on the coefficient ma- 
trix A, namely that all of its N x N maximal minors be non-negative. Such matrices are 
called totally non-negative matrices, they represent the totally non-negative Grassmannian 
Gr + (Af, M) which forms a closed subset in Gr(A r , M), and which identifies the space of 
all line-soliton solutions of the KP equation. 

Some simple examples of r-functions corresponding to some of the solutions men- 
tioned in Section 1 are listed below. 

One-soliton solution: In this case, one chooses a single function in © of the form 
fi = e® 1 + e 6 * 2 . The coefficient matrix is simply A = (1 1), and r = fi. Then equation 
(0) yields the one-soliton solution given by Q in Section 1 with # = 6*02 — #01 • Note 
that since ki < fo, r ~ e 6 * 1 3> e° 2 and u(x, y, t) ~ as x — > —00. A similar argument 
shows that u is also exponentially small as x — > 00. In this case, one finds from the exact 
expression in © that the solution is localized in the xy-plane along the line 61 = 62 (for 
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fixed t) where the two exponentials in the r-function are precisely in balance (see Figure 
1(a)). It turns out that even for general line-solitons, the solution u(x,y,t) is localized 
along certain lines in the xy-plane where exactly two exponential terms in the r-function 
of equation (0) are in balance, and they dominate over all other exponential terms. This 
principle of dominant balance was implemented in Refs 12111 to identify the 



asymptotic line soli tons [i, j] associated with each line-soli ton solution of the KP equation. 

Miles Y-shape solution: The r-function for this solution is a sum of three exponential 
functions, i.e., r = f\ = e 01 + e 02 + e 03 with the coefficient matrix A = (1 1 1). 
Applying the dominant balance principle mentioned above, it is possible here to determine 
the dominant exponentials and analyze the structure of the solution in the xy-plane. One 
finds that the solution u(x,y,t) is localized along three lines corresponding to the line 
soliton [1,3] for y 3> 1, and line soli tons [1,2], [2, 3] for y <C 1 as illustrated in Figure 
1(b). For i, j G {1,2,3}, each line soliton i < j is locally of the form of the 
one-soliton solution (0 with distinct parameters ki, kj. As mentioned in Section 1, the 
Miles solution represents a resonant solution of three line solitons. The resonant condition 
among those three line-solitons is given by 



43 



ki2+k23, 



where = |(fcj — ki,kj — kf) and = |(&| — kf) are respectively, the wave vector 
and frequency of the line soliton satisfying the soliton dispersion relation given in 
Section 1. 

O-type 2-soliton solution: Here the r-function is a Wronskian of two linearly indepen- 
dent solutions of (f5]). These solutions are obtained from ([6]) by choosing N = 2, M = 4, 
and the coefficient matrix 



.4 



110 
11 



such that fi = e 01 + e 02 and fi = e 03 + e 0i . Note that A is a totally non-negative matrix 
whose 2x2 maximal minors are: A{\,2) = A(3,4) = 0, and .4(1,3) = A(l,4) = 
A(2, 3) = .4(2, 4) = 1. From ([7]), the resulting r-function is given by 

r = (fc3-A;i) e e ( 1 ' 3 )+(A;4-fci) e e ( 1 ' 4 )+(fc3-fc2)e e(2 ' 3) +(A : 4-A ; 2)e e ( 2 ' 4 ) , B(i,j) := Oi+Oj 

Figure 1(c) illustrates the regions in the xy-plane where each of the four exponential terms 
in the above r-function are dominant. The line solitons corresponding to the solution 
u = 2(Iiit) xx are [1,2] and [3,4] for \y\ — > oo; these are localized along the directions 
where a pair of dominant exponential terms are in balance. 



2.2. The Grammian form of the T-function 

The Grammian construction of the line-soliton solutions will be described next. This form 
of the solution arises from the so called binary Darboux transformation of the KP equation 



(see e.g., [24]). In this case, the r-function can be expressed as a a determinant of an NxN 



matrix as follows: 



det(fi) , ftij = J figj dx-{fig jx - f ix gj) dy , i,j = 1, 2, . . . , N , (8) 
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where {fi]f =l solve © and {gj}j =1 solve the (formal) adjoint system of ©, namely, 

9y 9xx i 9t 9xxx • 

It is easily verified from (HJ) and its adjoint that dQij is an exact differential, hence the line 
integral above can be evaluated along any suitable curve T in the xy-plane such that the 
integral converges. Like the Wronskian form, it can be shown that the r-function in © 
yields a solution of the KP equation via (0) for any choice of the sets of functions {fi}jL 1 
and {gj}f =1 & For the line-soliton solutions, these are chosen as linear combinations of 
exponentials (cf. ©), 

r N 

fi = Y, h ^ n \ 9 j = ^2b j ne- (t>M , i,j = l,...,N, (9) 

n=l n=l 

where cj)(k) := kx + k 2 y + k 3 t + 4>o(k), and all the parameters pi, . . . ,p r , qi,...,qN 
are distinct and real. Note that N 7^ r, in general. Thus, B := (b mn ) is an N x r matrix 
whereas B := b mn is an invertible, N x N matrix. Substituting © into the matrix ele- 
ments Qij in ([8]), and evaluating the line integral along a path F from (—coo, y) to (x, y) 
parallel to the x-axis, yields 

r ' _^ ^ e <t>(p m )-<t>(<in) 
£lij(x, y, t) = dj + / fiQj dx = Cij + 2^ 2^ bimbjn , (10) 

jT 1 1 Pm qn 

m=l n=l 

where a = sgn(p m — q n ) = ±1 and are arbitrary constants. Choosing the constants 
Cij = 5ij, where <5y is the Kronecker symbol (5ij = 0, i 7^ j, ^ = 1), ensures that the 
matrix 17 is of rank N so that r^O. Using (fTOl ). the r-function in ([8]) can be expressed in 
the form 

r = det(J + BFB T ) = det(J + CF) , (1 1) 

where C = B T B is an iV x r matrix of constant coefficients and F is an r x N matrix 
whose entries are given by 

Fmn = , m = l,2, ...,r, n = 1,2, . . . ,N. 

Pm qn 

The Grammian form given by the first equality in (fTTT > is not unique as it is possible to 
obtain the same r-function from a different choice of the matrices B and B. In particular, 
B need not be a square matrix although there is always a canonical choice of an N x N 
matrix B that leads to the second equality in (TTTb . We will not discuss the details of the 
gauge freedom underlying the choice of B and B in this article. 

Equation (TTTb is the canonical Grammian form for the line-soliton solutions, special 
cases of which arise from several direct methods of constructing solutions including the 
Hirota method |@], direct linearization 10] and dressing techniques For example, the 
O-type iV-soliton solution is obtained by setting r = N and C = I, the latter being the 
N x N identity matrix. In this case the resulting r-function r = det(I + F) is positive for 
all x, y, t if the parameters are ordered as < pn-i < ■ ■ ■ < Pi < <Zi < qi < ■ ■ ■ < qN- 
Then the corresponding line-soliton solution is non-singular. More general line-soliton 
solutions can also be generated from the Grammian form (TTTb by making appropriate 
choices for the matrix C. However, in comparison with the Wronskian form, here it is 
less clear how to impose regularity conditions on the obtained solution. Recall that in 
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the Wronskian construction the ordering k\ < k 2 < • • • < &m of the parameters and 
the totally non-negative coefficient matrix A guarantee that the KP line-soliton solutions 
are regular. But no such clear prescription to obtain regular solutions is known for the 
Grammian case. This issue can be resolved by establishing an equivalence between the 
Wronskian and the Grammian forms for the KP T-function. This will be done below. 
More specifically, it will be shown by explicit construction that there is a one-to-one 
correspondence between the formulas |7]) and (fTTT) . 



2.3. Equivalence of the Wronskian and Grammian forms of the T-function 

It is convenient to first start with the Wronskian form and derive the r-function formula 
CED from that. It follows from either (0 or (@) that the r-function can be written as: 
t = det(AEK), where E = diag(e ei , e 02 , . . . , e 0M ) is an M x M diagonal matrix 
and K := (Kij) = (kj' 1 ) is an M X N Vandermonde matrix. In the following, the 
N x M coefficient matrix A will be chosen in RREF such that it can be represented as 
A = (I , J)P with I and J denoting respectively, the N x N and N x (M — N) sub- 
matrices of the pivot and non-pivot columns of A, and P being the M x M permutation 
matrix which shuffles those columns to form A. For example, the coefficient matrix for 
the O-type 2-soliton r-function discussed above can be represented as 

. (\ 1 OCA /10 10\ „ 

A ~ iooiij " VoioiJ P23 ' 

where P23 is the 4x4 matrix permuting the second and the third columns. Using the 
above form of A, the determinant for the r-function becomes 

r = \AEK\ = \(I,J)(PEP- 1 )PK 

where E\ and E2 are respectively, N x N and (M — N) x (M — N) block diagonal 
matrices whose elements are permutations of the set {e e% , e® 2 , . . . , e eM }. Similarly, K\ 
and K2 are respectively, N x N and (M — N) x N matrices obtained by permuting the 
rows of the Vandermonde matrix K by P. It should be clear from above that in effect, the 
M x M matrix P induces a permutation ir of the ordered set {ki, k%, ■ ■ ■ , &m} which can 
be expressed as 

TT({k 1 ,k 2 ,... : k M }) = {gi,<?2, ...,qN,Pi,P2, ...,Pm-n}, (12) 

after renaming the elements of the permuted set. Accordingly, the matrices 
Ei , E2 , K\ , if 2 are redefined as 

Eh = diag(e^ qi \e^ q2 \. . . , e^) , E 2 = diag(e^ Pl \ e^ P2 \ e ^ p "-^) , 
Ki := {K Xij )» j=1 = , K 2 := (K 2ij ) = {^) , 

for i = 1, . . . , M — N, j = 1, . . . , N, and where (j)(k) is defined below ©. Since the 
parameters q\, . . . , are distinct, the Vandermonde matrix K\ is invertible. Then the 
above determinant expression for r can be further manipulated as 

r = \E x K x \r , ?=\I+ JE 2 K 2 K^E^ l \ , 

where r and r differ by an exponential factor linear in x,y,t, hence both generate the 
same KP solution via ©. To show that r is indeed the Grammian form of the r-function, 



(I, J) 



Ei 
E 2 



Ki 
K 2 



\IEiKi+JE 2 K 2 \ 
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one employs the following matrix identity which is derived in the Appendix, 

K 2 K~ X = D 2X D^ , (13) 
where the matrices D\,D 2 and the Cauchy matrix X are defined as follows: 

£>i = diag^ [[ (Qi-Qj)J , D 2 = diag(^ [[(pi - qj)j , 

X = (x)ij = — , i = 1, • • • , M - N, j = 1 . . . , N . 
Pi ~ Qj 

Substitution of ( [TBI into the expression for r above, yields 

t = \I + JD 2 E 2X E^D^\ = \I+ {D^ l JD 2 ){E 2X E^ l )\ . 

Finally, by setting r = M — N and the N x M — N matrix C = D^ l JD 2 , the Grammian 
form of the r-function in (fTTT) is readily recovered from the expression of r above. 

It is relatively straightforward to reverse the steps described above to obtain the Wron- 
skian data, i.e., the coefficient matrix A and the ordered set of parameters k\ < k 2 < 
. . . < ku from the Grammian data, which consists of the N x M — N matrix C and 
the unordered set of parameters {qi , . . . gjy, pi , . . . pm-n}- The key step is to recover the 
permutation tt in (fT2l ) by simply ordering the set of Grammian parameters, i.e., 

TT~ 1 ({q 1 ,...q N ,p 1 ,...p M ~N}) = {h,k 2 , . . . ,k M } ■ 



This in turn, provides the M x M permutation matrix P from tt, and the matrix A is 
then constructed by the formula A = (I , J)P with J = DiCD^ 1 , thus establishing the 
Grammian-Wronskian equivalence for the r-function for the KP line-soliton solutions. 
Such equivalence was also derived in Ref. J25ll for special iV-soliton solutions using sim- 
ilar techniques. 

An alternative approach to obtain the Wronskian form from the Grammian in ([8]) that 
is applicable when at least one of the functions fi or gj is an exponential, was discussed 
in 12611 . In this case, one can for instance, choose the functions {/j}^ as N linearly 
independent solutions of ©, and {gj := e~ 6j }f = i with 9j = kjx + k 2 -y + kjt and kj > 0. 
Using these forms of fi and gj in ([8]), and choosing T to be the path from (— 00, y) to (x, y) 
parallel to the x-axis, the integral for the matrix element f^- can be expressed as an infinite 
series in inverse powers of the parameters kj by repeated integration by parts, provided 

(n) 

that /> gj —> 0, n > as x — > —00. Then by taking the limit kj — > 00, j = 1, . . . , N, 
one can recover the Wronskian form ((U) from the determinant det(fi). 

An immediate consequence of the above equivalence result is that it is now possible 
to formulate the regularity condition for the KP line-solitons in a precise fashion from 
the Grammian form (fTTT) of the r-function. Given the matrices C and P, a necessary 
and sufficient condition that the KP solution (0 is non-singular is given by the fact that 
A = (7, J)P be a totally non-negative matrix. Of course, this is the same condition on 
the coefficient matrix A in the Wronskian form. If in addition, A satisfies the following 
irreducibility conditions: 

(i) each column of A contains at least one nonzero element, 

(ii) each row of A in RREF contains at least one nonzero element other than the pivot 
(first non-zero entry), 

then the number of line soli tons as y — > ±00 in the solution generated by the r-function 
in ([7]) is determined by the size of the matrix A. Namely, one has A r _ = M — N and 
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N + = N, where N± denote the number of line solitons as y — > ±00 II20U2 11 12311 . Further 
analysis of the totally non-negative matrix A leads to the precise identification of the line 
solitons and a comprehensive classification scheme for all line-soliton solutions of the 
KP equation. The latter problem is related to the classification of the non-negative cells 
Gr + (iV, M) of the Grassmannian Gr(iV, M) in terms of certain types of permutations 
called derangements J2lll . 

It is worth noting that the equivalence between the Grammian and Wronskian forms 
of the line-soliton r-function is not unique to the KP equation but applies also to other 
equations whose r-functions are represented in both forms. Examples include several 
2 + 1-dimensional differential-difference equations such as the 2d-Toda, 2d-Volterra and 
the differential-difference KP equation as well as difference-difference equations such as 
the fully discrete 2d-Toda and the discrete KP (Hirota-Miwa) equation. These will be 
discussed in future works. 



3. O- and P-type 2-solitons 



In this section the soliton interactions associated with certain types of 2-soliton solutions 
of the KP equation will be investigated. These solutions are characterized by a pair of line 
solitons as \y\ — > 00. Thus, iV_ = N + = 2, which implies that N = 2 and M = 4. The 
Wronskian form of the r-function for these solutions is obtained from ^} in terms of the 
ordered set k\ < k 2 < &3 < &4 of distinct real parameters and a 2 x 4 totally non-negative 
matrix A as, 

r(x,y,t)= (ks-k r )A(r,s)e e ^ , (14) 

l<r<s<4 



where A(r, s) denotes the 2x2 minors of the matrix A. Taking into account the irre- 
ducibility property mentioned in Section 2, there are exactly seven distinct types of totally 
non-negative 2x4 matrices, which lead to seven different types of 2-soliton solutions 



of the KP equation 112 11 123H . In contrast, its (l+l)-dimensional version namely, the KdV 
equation has only one kind of 2-soliton solution. This is indicative of a much richer so- 
lution space for the (2+l)-dimensional integrable equations than the (l+l)-dimensional 
ones. 

The particular 2-soliton solutions whose nonlinear interaction properties are studied 
here, are the O-type 2-soliton which was mentioned in section 2, and another one called 
the P-type ("P" for physical), which fits better with the physical description of oblique 
interactions of shallow water waves for which the KP equation is a good approximation. 
In both of these cases, a pair of line solitons interact to form a X-shape in the xy-plane but 
the details of the interactions are entirely different. The corresponding coefficient matrices 
in RREF are given by 

/110 (A . A 0-1 

A ° ~ [0 1 1) ' Ap ~ [p 1 1 

which, from (fT4b . lead to the r-functions 

r = (fc 3 - kt)e e ^ + (fcg - k 2 )e 9 ^ + (k A - h)e ^ + (k 4 - k 2 )e 9 ^ , 
Tp = (fe 2 - ky)e e ^ + (k 3 - kx)e 9 ^ + (k A - k 2 )e ^ + (fe 4 - k 3 )e e ^ . 



Applying the principle of dominant balance mentioned in the examples given in Section 2, 
it is possible to identify the line solitons in the two cases. For the O-type the line solitons 
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Figure 2. An O-type 2-soliton. The interaction peak corresponds to an absolute maximum. 



are [1, 2] and [3, 4] as y — > ±oo, whereas for the P-type, these are [1, 4] and [2, 3]. Recall 
that locally, a line soliton [i, j] is given by ([2]) with parameters hi and kj. 

For both O- and P-type solitons, the exact solution can be computed as 



u = 2(lnr) 



a sech 2 <f>x + /3 sech 2 <fi 2 + j seen 2 0x sech 2 <f> 2 
1 + C tanh^>itanh^>2 



with a = A^l - C 2 ), (3 = A 2 {1 - C 2 ), and 7 = (Ax + A 2 )C 2 + 2CjA~x~A 2 . The 
expressions for the soliton amplitudes A\,A%, soliton phases 4>i,4> 2 , and the constant C 
are different for the O- and P-type 2-solitons. These are given below. 



3.1. O-type interaction 

For the O-type, the amplitudes of the line solitons [1, 2] and [3, 4] are 

Ax = a 12 = \{k 2 - ki) 2 , A 2 = q 34 = \{k A - k 3 ) 2 , 

the soliton phases are expressed in terms of Oj = kjx + k 2 y + kjt + 6oj, j = 1, . . . , 4, as 



h = ^(0 2 -0i)+Tln 



(fc 4 ~ k 2 ){h ~ k 2 ) 
(h - h)(k 3 - ki) 



(fe 4 - k 2 )(k A - fci) " 
(h ~ k 2 )(k 3 - kx) 



and the constant C is given in terms of the parameters kx , . . . k^ by 



C 



o 



1 + VA 



A 



o 



o 



(^3 ~ fe2)(^4 - fcj) 

(k 3 - k 1 )(k 4: - k 2 ) 



(k 2 - fci)(fc 4 - k 3 ) 
(k 3 - kx){k/x - k 2 ) 



< 1 



such that < C < 1. An O-type soliton is illustrated in Figure 2. The function u(x, y, t) 
has an absolute maximum at cf)x-<j) 2 - 0, and the maximum value at this interaction point 
is given by (see l23ll ) 



Ax + A 2 + 2Cy / A 1 A 2 . 



(15) 



Since C > 0, the interaction peak u max is always greater than the sum of the individual 
line soliton amplitudes. In fact, one can easily verify that 



Ax + A 2 < u max < ( V^i + \[A 2 f . 
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Furthermore, u max depends non-linearly on the incidence angle between the two line 
solitons. In order to investigate this behavior, it is convenient to choose the parameters 
fei, . . . , &4 such that k\ = —k^ and k2 = —k% = —k (see Figure 3(a)). Then the two 
line solitons are placed symmetrically about the y-axis, and are of equal amplitude, i.e., 
A\ = A2 = A. Denoting the angle between the [3, 4]-soliton and the y-axis by ip as shown 
in Figure 3(b), it follows that tan ip = (k^ + k^) = y/2A + 2k. Since k > 0, the angle ip 
is always greater than the critical angle 



*J\: 



tan 



f 2A) 



which depends on the amplitude of the line solitons. Setting k\ 
expression for C above, it can be deduced from (PT51) that 



-ki, ko 



-fe in the 



4A 



i + VT 



tan ip c 
tan ip 



< 1 



Therefore, u max is a decreasing function of the angle ip, and as ip — > ip c (i.e., k — > 0), the 
peak interaction amplitude n max — > AA. Note that in the limit k = 0, the parameters /c2 = 
&3, which implies that the r-function to contains only three instead of four exponential 
terms. The resulting solution u(x, y, t) of the KP equation is a Miles' Y-shape solution 
corresponding to the confluence of three line solitons interacting resonantly 111 ill (see also 



123Q ). The dependence of n max with the angle of interaction is shown in Figure 4 for 
soliton amplitudes A = 1 and A = 2. Note that in Figure 4(a), u max — ► 4 as ip approaches 
the critical angle ip c 54.7°, while in Figure 4(b), u max — > 8 as ip — ► ipc ~ 63.4°. 
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Figure 4. Plots of u max vs ip for O-type 2-solitons. The amplitudes of the line solitons [1, 2] and [3,4] are fixed and 
equal, i.e., Ai = A 2 = A: (a) A = 1, ip c « 54.7°; (b) A = 2, ip c a 63.4°. 
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(a) 

Figure 5. (a) A P-type 2-soliton. The interaction between the line solitons is characterized by a saddle point; (b) Contour 
plot of a P-type 2-soliton with fci = — 3, &2 = —0.5, k% = 2, &4 = 3. The line soliton [1, 4] is along the y-axis, the angle 
between the line solitons [1, 4] and [2, 3] is ip. 

3.2. P-Type interaction 

For the P-type, the soliton amplitudes 

Ai = au = \(fa - fa) 2 , A 2 = a 2 3 = \(fa - fa) 2 ) 

are always unequal, and A\ > A 2 . The soliton phases are given by 



(fa - fa) (fa - fci) 
(fa - fa)(fa - fa) 



2 y J 4 
and the constant C is negative since 
1 - v'Ap 



l -(e z -9 2 )+-\n 



(fa - fa)(fa - fa) ' 
(fa ~ fa)(fa -fa). 



c 



1 + 



= (fa - fa)(fa - fa) 
(fa - fa)(fa - fa) 



A P-type 2-soliton solution is illustrated in Figure 5. In contrast to the O-type, the P-type 
2-soliton solution u(x, y, t) corresponds to saddle at the interaction point 4>\ = 4> 2 = 0. 
The value of the solution at the interaction point is 



A x + A 2 + 2C y /A 1 A 2 ■ 



(16) 



Since C < for the P-type solution, u- mt is always less than the sum of the amplitudes of 
the line solitons [1, 4] and [2, 3]. A more precise bound for u- m t is obtained as follows J23ll : 
From the above expression for C, one can calculate 



C 



1 



7(l-x)(l-y) 



xy 



1 - (x + y) 



where x = and y 



l + V^P y/(l-x)(l-y) + ^xy W(l-x)(l-y) + Jxy)2 
k^zzhs si nce q < x, y < 1, the following inequalities 



xy < Mx + y) 



y/(l-x)(l-y) <l-Ux + y) 



hold, and from these it can be easily deduced that 1/(1 — x)(l — y) + ^/xy < 1. The 
latter inequality leads to 



1 - ^Ap / . \ fa -fa 
■ > 1 — (x + y) = 



fa — fa 
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Thus, -1 < C < -y / A 2 /Ai, and from £l6]> it follows that 



Next we investigate the dependence of the amplitude u- mt at the interaction point on the 
angle between the two line solitons [1, 4] and [2, 3] with fixed amplitudes Ai and A 2 . For 
simplicity, we take the [1, 4] -soli ton along the y-axis by setting k\ + = 0, and denote 
by ip the angle between the solitons as in Figure 5(b). Then 



tan -0 = k 2 + k 3 = 2{k 2 - h) - {y/2A~i - \f2A 2 



after using k\ + k^ 
which implies that < k 2 



0. Similarly, one calculates 
< \/2A~i 



h) + (k 4 - h) = y/2Al- ^/2A 2 



k\ < \/2A\ — \J2A 2 . Hence for fixed amplitudes A\ and 
A 2 , the angle ip between the two line solitons satisfies 



-ip c < ip < 1p c 



ip c = tan"^ yj2M - y/2A 2 ) 



The limiting cases ip = ±ip c correspond to k\ = k 2 or fc 3 = fc 4 . Either case leads to 
a degeneration of the P-type 2-soliton to a Y-shape solution similar to O-type 2-soliton 
situation. In terms of A±, A 2 and the angle ip, the quantity Ap can be expressed as 



A, 



(k 2 - h)(k 4 - k 3 ) 



tan 2 ip c — tan 2 ip 



and then from ( fT6l ) one obtains 



k 2 ) {y/2A^ + y/2A 2 ') 2 - tan 2 ip 



Uint 



tan 2 ip + \J (tan 2 ip c - tan 2 ip)((y^2Al + ^/2A 2 ) 2 - tan 2 ip) 



It follows from the above expression that u- mt is an even function of the angle ip between 
the line solitons [1, 4] and [2, 3] for \ip\ < ip c . When ip = 0, i.e., when the two line solitons 
are both parallel to the y-axis, u- mt reaches its maximum value of A\ — A 2 which is the 
difference between the amplitudes of the two line solitons. As ip — > ±ip c , the interaction 
amplitude approaches its lower bound, i.e., u- mt — > ^ tan 2 ip c = (\/A~i — \/A 2 ) 2 . Figure 
6 below illustrates the plots of u- mt as a function of the angle ip for fixed values of soliton 
amplitudes A\ and A 2 . 
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Figure 6. Plots of u- lnt vs interaction angle ip for P-type 2-solitons with fixed amplitudes Ai, A2 of the line solitons 
[1,4] and [2,3]: (a) Ai = 2, A\ = 1 and ip c 30.4°. The maximum value of u- mt is 1 at ip = 0, and u; nt — » 
(\/Ai — V 'A2) 2 ~ 0.17 as ip — » ±i/> c ; (b) A\ = 4, A\ = 1 and ip c 54.7°. The maximum interaction amplitude 
(tiini)moi = 3 at ip = 0, and u int — > 1 as ip — > ±ip c . 
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4. Conclusion 

In this article we have considered a special class of non-singular solutions of the KP equa- 
tion referred to as the line-solitons which decay exponentially as x, y — > ±00 except along 
certain directions in the xy-plane. Such solutions exhibit a variety of time-dependent spa- 
tial patterns due to resonant soliton interactions (Y-shape) as well as non-resonant interac- 
tions (X-shape). The exact analytic form of these solutions can be derived from the asso- 
ciated r-functions which are expressible either as Wronskians or as Gram determinants. It 
is remarkable that the KP equation possesses such a rich structure of line-soliton solutions 
generated from a simple form of the r-function. It turns out that the solution manifold of 
the line-solitons is parametrized by a discrete set of real distinct parameters and the space 
of totally non-negative matrices. This characterization is clear from the Wronskian form 
of the KP r-function but not so transparent from its Grammian form. It is perhaps due 
to this difficulty in imposing appropriate regularity conditions that only a handful of line- 
soliton solutions were explicitly known via the direct algebraic methods which used the 
Grammian form of the r-function. This issue has been resolved in this article where a 
one-to-one correspondence between the two forms of the line-soliton r-function has been 
established. Consequently, it is now possible to derive non-singular line-soliton solutions 
using the Grammian form of the r-function as well. 

Another problem discussed in this article is the nonlinear soliton interactions for certain 
types of 2-soliton solutions. In particular, the amplitude of the nonlinear interaction has 
been explicitly calculated from the exact analytic expression for each of the 2-soliton 
solutions considered in this paper. Moreover, the dependence of the interaction amplitude 
on the angle between the line solitons has been investigated by keeping the line soliton 
amplitudes fixed. One possible physical significance of such results lies in the study of 
oblique nonlinear interactions of weakly 2-dimensional solitary waves in shallow water. 
The physical mechanism generating large amplitude waves of extreme elevations from the 
interaction of two (or more) smaller amplitude solitary waves in shallow water constitutes 
an important open problem. It is believed that in appropriate parameter regimes the line- 
solitons of the KP equation can serve as a reasonably good test-bed for the description 
and analysis of nonlinear solitary wave dynamics. As a qualitative evidence, one may 
consider the simple example of the O-type 2-soliton interaction where the interaction peak 
amplitude may reach as high as four times the amplitude of the individual line solitons. 

For both O- and P-type solutions, the range of this interaction angle is found to be lim- 
ited by a critical value ip c which depends on the amplitudes of the line solitons. It is how- 
ever important to note that the theory of KP line-solitons can still be applied to study wave 
interaction where the angle between the incident waves is outside of the prescribed range 
for the O- or the P-type soliton solutions. In such cases, the wave dynamics is g overned 
by other types of KP 2-soliton solutions which have recently been uncovered \ lJi 2(1 21]. 



The study of the interaction properties of these newly found solutions is a topic of future 
investigation. 
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Appendix A. 

Here we provide an elementary derivation of the matrix identity ( fT3T ) involving the Van- 
dermonde matrices K\ and A'2. Recall that K\ is an N x N matrix with elements 
{K\)ij = q{~ 1 , and K2 is an /V x (M — N) matrix with elements (-^2)13 = P^ 1 where 
{qi, <?2, • • • , qN,Pi,P2, ■ ■ ■ ,Pm~n} is a set of distinct elements. 

Consider N monic polynomials P%(x), P2(x), . . . , Pn{x) such that each polynomial is 
of degree N — 1 and has distinct real roots which consist of N — 1 elements from the set 
{qi, <?2, • • • , <?tv}- These polynomials can be expressed as 

N N 

Pj(x)-= II (x-Qk) = ^Pn J x n - 1 , j = l,2,...N, (Al) 

k=l,k^j n=l 

where P n j is the coefficient of x n ~ l in the polynomial Pj (x), and P^j = 1. An immediate 
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consequence of (IA 1 b is that the polynomials satisfy 

N 



p j(li) = Pnjli' 1 = ( K i p h = {Dihtij , h k = 1, 2, . . . , N , (A2) 

71=1 

where P = (P n j)nj-i * s tne matrix of polynomial coefficients, D\ is the diagonal matrix 
defined below ( fT3l ), and <5y is the Kronecker symbol. Equation (IA2b gives the matrix 
equation K\P = D\ which implies that 

P = K^ 1 D 1 (A3) 

since K\ is invertible. Next, evaluating the polynomials in (IA 1 b at x = pi, yields 

N 



N N 



no* 



PM) = Y, p njpr 1 = n fe-^) = ra=1 p ._ q . , 

n=l k=l,k^j P% tfjf 

for i = 1, . . . , M — N, j = 1, . . . , N. These can be represented by the matrix equation 

K 2 P = D 2X , (A4) 

where the diagonal matrix D 2 and the Cauchy matrix x are a l so defined below (fT3l . 
Combining (IA3b and (IA4b gives the desired identity (fl"3T ). 



